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Kafri, Mukamel, and Peliti Reply: The Comment 
of Hanke and Metzler [jij questions the validity of the 
analysis presented in Q to DNA chains of finite length as 
used in experiments. Their argument is that for the anal- 
ysis to be valid “each of the three segments going out of a 
vertex must be much longer than the persistence length 
l p of this segment”. By using the persistence lengths 
£ P (L) ~ 40Å for a single strand and £ P (H) ~ 500 Å for 
a double helix (bound segment) they arrive at the con- 
clusion that in order to observe the asymptotic behavior 
found in j^j] one needs chains which are far longer than 
those studies experimentally. 

This assertion constitutes a misunderstanding of the 
analysis given in Q . In this analysis one considers a loop 
interacting with the rest of the chain and not just with the 
vicinal double helices. Thus, in [Q each of the two lines 
attached to the loop is in faet composed of an alternating 
sequence of bound segments and denaturated loops. It 
may be viewed as a stick and joint structure, whereby 
adj acent double helices (which may be considered as rigid 
rods as they are shorter than £ P (H)) are loosely attached 
to each other via an open loop (see Fig. 1). Thus, the 
“rest of the chain” as considered in reference [2j is in 
faet a random rod structure with a persistence length 
given by the length of the rod connecting two loops. This 
length can be easily estimated from the analysis of Poland 
and Scheraga and its generalization as given in |l|. It 
can be shown that the probability distribution of rods 
of length k at the transition is P(k) ~ ( w/s) k . Here 
w = exp(/3M-Eb) is the Boltzmann weight of a bound 
pair with energy —Eq at the melting temperature, 1//3m, 
and s is a non-universal geometrical factor. This is an 
exponentially decaying distribution with a typical length 
of order £ ~ 1/| ln(w/s)|. This length has nothing to do 
with the persistence length of the bounded segments and 
in faet it is far shorter (of the order of a small number of 
base pairs). An estimate of this length can be obtained 
using Eq. 4 of [Q at the transition 



With the estimated value of c = 2.115 (in three dimen- 
sions) this yields £ cs 1.07 in units of £ P (L). Note that 
this length is non-universal and it depends on the de- 
tails of the model, such as the statistical weight of small 
loops, the stacking interactions and other details. How- 
ever, these features are not expected to change this length 
considerably. Therefore, the criticism expressed in the 
Comment is not valid and the analysis in the Letter is 


applicable to finite chains of the length studied experi- 
mentally. In reference || the detailed structure of the 
“rest of the chain” is ignored when the self-avoiding in- 
ter action is considered. However, as ar gued in Q this 
assumption is reasonable and should yield a good esti- 
mate of the effeet of self-avoiding interactions. 



FIG. 1: Schematic representation of a microscopic configura- 
tion of the DNA molecule. Here the flexibility of the “ rest of 
the chain ” emerges from the presence of loops which connect 
short but more rigid bound segments. 


A strong support for this picture is provided by nu- 
merical simulations of Causo et. al. || and Carlon et. 
al. [Q of chains of finite length in which the self avoiding 
interactions are fully taken into account. In both studies 
the melting transition was found to be first order, com- 
patible with c > 2 as found in (Q . In Q the exponent c is 
evaluated direetly for chains of length 50-200 monomers 
(each monomer corresponds to £ P {L) namely to about 
8 base pairs). Moreover, the inclusion of the different 
stiffnesses of the bound and unbound segments (in simu- 
lations persistence length ratio, £ p (H)/£ p (L) : was taken 
to vary between 1 and 10) did not change the estimate 
significantly. The faet that the observed c is consistent 
with that obtained by reference |J strongly supports the 
basic assumptions behind this analysis. It also indicates 
that the scaling behavior predicted in Q] may be observed 
in finite chains of lengths between hundreds to thousands 
base pairs, which is the experimentally relevant range. 
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